Multifractality at the spin quantum Hall transition 
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Statistical properties of critical wave functions at the spin quantum Hall transition are studied 
both numerically and analytically (via mapping onto the classical percolation). It is shown that 
the index rj characterizing the decay of wave function correlations is equal to 1/4, at variance with 
the r~^^^ decay of the diffusion propagator. The multifractality spectra of eigenfunctions and of 
two-point conductances are found to be close-to-parabolic, A, ~ g(l — q)/S and Xq ~ g(3 — g)/4. 
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Disordered two-dimensional electron systems show re- 
markably rich physics, which is governed by quantum 
interference effects and depends on the symmetry class 
to which the system belongs. Recently, unconventional 
symmetry classes Q , which can be realized in d-wave su- 
perconductors, have attracted considerable interest. Par- 
ticular attention was paid to the class C (broken time- 
reversal invariance) , where a transition between localized 
phases with quantized values of the spin Hall conduc- 
tance takes place . A network model describing this 
spin quantum Hall (SQH) transition was constructed in 

, and critical exponents for the scaling of the localiza- 
tion length were determined numerically. In Q| a map- 
ping onto a supersymmetric spin chain was performed, 
providing an alternative method for the numerical study 
of the critical behavior. Remarkably, some exact analyti- 
cal results for this model have been obtained by mapping 
onto the classical percolation problem . In particular, 
it was found that the density of states (DOS) is critical 
and scales as p(e) ~ e^/^, while the diffusion propagator 
H(r,r') = (Gi<;(r, r')GA(r', r)) and the average two-point 
conductance (g{r,r')) fall off as |r — r'|~^/^ at criticality. 

The aim of the present paper is to study the statistical 
properties of wave functions at the SQH critical point. 
Multifractality of wave functions ip{r) is known to be a 
hallmark of the localization transition. It has been ex- 
tensively studied in the context of conventional Ander- 
son and quantum Hall (QH) transitions with non-critical 
DOS (see e.g. ^ and references therein), and we remind 
the reader of some basic results. Multifractality is char- 
acterized by a set of exponents Tq = d{q — 1) + Aq (d is 
the spatial dimensionality) describing the scaling of the 
moments of |7/'^(r)| with the system size, (IV'(r)P^) oc 
Ij-d-Tg Anomalous dimensions Aq distinguish a critical 
point from the metallic phase and determine the scale 
dependence of wave function correlations. Among them, 
A2 = —77 plays the most prominent role, governing the 
spatial correlations of the "intensity" , 



functions and the diffusion propagator possess the same 
scaling properties, 

L2''(|V','(r)^|(r')|), L^'mv)r,{v)rd^')U^')), 
p-^H(r,r';^)^(|r-r'|/i.)-^ (2) 



where w — ei — ej, ~ (pw)"-'^/'' and |r — r'| < L^. The 
moments of the two-point conductance show a power-law 
scaling as well [|^,^, {g''{r,r')) oc |r — r'|^'''-« with another 
set of exponents Xq, which are related to A^ [p|,fol. 
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Correlations of two different (but close in energy) eigen- 



In two dimensions the multifractal spectra Aq and Xq 
play a key role in the identification of the conformal field 
theory governing the critical point, which led to growing 
interest in the eigenfunction statistics at the QH transi- 
tion 

Applying naively these results to the SQH transition, 
one would conclude that the r^^/^ scaling of the diffusion 
propagator found in implies rj ~ 1/2. However, a more 
careful examination of the findings of Ref. Q gives us a 
warning: Gruzberg et al. find that H and {g) share the 
same r~^/^ scaHng, in contrast to Eqs. (||), (||) predicting 
different exponents (77 = — A2 and 2A]^/2j respectively). 
This suggests that one should be cautious when trying to 
apply the above relations between critical exponents ob- 
tained for systems with a non-critical DOS to those with 
a critical one (like the SQH transition). Our analysis 
below shows that such a caution is well justified. 

After this introduction, we turn to our study, which 
combines numerical and analytical methods. Our main 
aim is to calculate exponents governing the scaling of 
dV'l^''), {g''), and H, and to understand relations between 
them. For computer simulations we used the SU(2) net- 
work model described in Each realization of the 
network is characterized by a unitary evolution operator 
lA. Diagonalizing lA using advanced sparse matrix pack- 
ages yields eigenfunctions ijji and eigenvalues e~"', 
where i = 1,2,...L^. In Fig. |l] we display the DOS 
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for different system sizes L. It is seen that after a proper 
rescaling all data collapse onto a single curve, as expected 
at criticality. At e 3> 5 (where 5 '-^ L^^/^ is the level spac- 
ing at e = 0) the DOS scales as p(e) ^ e^/^, in agreement 
with the analytical prediction On the other hand, at 
e ~ 5 one observes an oscillatory structure qualitatively 
analogous to the behavior found in the random matrix 
theory for the class C 
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FIG. 1. Scaling plot of the density of states for system sizes 
L — 16(o), 32(Q), 96(o). Dashed and dotted lines indicate 
power laws (dashed: e^^'^ , dotted: e^), S — 1/2-kL7^'^ denotes 
the level spacing. Inset: same data on a linear scale and the 
result from random matrix theory [Q] (solid curve). 

We concentrate now on the statistics of eigenfunctions 
with smallest energies, |e| ~ 5, for which the correlation 
length ^ e^^/'' Q is of the order of the system size. 
In Fig. 1^ (inset) we plot the cigenfunction autocorrela- 
tion function (|l|), as well as the correlation function of 
two eigenstates neighboring in energy (the second one in 
Eq. The results imply a power-law behavior oc r^"^ 
of both correlation functions at distances 1 <C r ^ L, 
with an index i] close to 1/4. On the other hand, our nu- 
merics confirms the value 1/2 of the exponent govern- 
ing the decay of 11 (Fig. |[ main panel). To understand 
the difference between the two exponents, we turn now 
to an analytical approach. 

Consider a correlation function of two wavefunctions, 

V{e',e;ei,e2) = i/'L(e)i/'ja(e)i/'«/3(e')V'j/3(e') 

xS{e,-e.)S{e2-e,)), (4) 

where e, e' are two edges of the network and a, f3 — 1,2 
are the SU(2) indices. Introducing the Green function 
G{e',e;z) = (e'|(l — zL{)~^\e), we express (jj) as 

V{e\ e; ei, €2) = (27r)-2(Tr[Gfl(e', e; - ^^(e', e; e'^^l 
x[Gfl(e,e';e'^^)-G^(e,e';e"=)]), (5) 



where Gr^a are retarded and advanced Green functions, 
Gii,A(e',e;e"i) = G(e', e; e*(<^i±^°)). We will calculate 
(^ at zero energy, £12 0, but finite level broadening, 
±iO — > ±47. The scaling behavior of the correlation func- 
tion (P at ei, £2 ~ e can then be obtained by substituting 
e for 7. We thus need to calculate 

P(e', e; 7) = (2^)-2(Tr[G(e', e; z) - G(e', e; z'^)] 

X [G(e,e';z)-G(e,e';z-i)]), (6) 

with a real z = e ' < 1. To do this, we make use of 
the mapping to the classical percolation, following the 
approach of Q . We give only a brief outline of the calcu- 
lation here; details will be published elsewhere ||l^. The 
Green functions in (H) are represented as sums over paths; 
the resulting expression is to be averaged over SU(2) ma- 
trices Uf associated with all network edges /. The crucial 
point is that for each edge / only paths visiting it or 2 
times are to be taken into account. In |^ this was proven 
for the average Green function (TrG(e, e; z)). The proof 
is based on the observation that (Uj) = Cq ■ 1, where 
Cq — for integer q ^ 0, ±2. We generalize the statement 
to products of two Green functions of the type entering 
as well as another two-point correlation function, 

D{e\ e; 7) = (2^)-2(Tr[G(e, e; z) ~ G(e, e; z'^)] 

X [G(e',e';z)-G(e',e';z-i)]), (7) 

corresponding to the (e)V'|(e')|) correlator, in the 
following way. Classifying the paths according to the 
number of times they return to a link /, we obtain ex- 
pressions of the type J2'^^^q^=i{TrUj' AUf B)xi^+i^ with 
A,Be SU(2) and x S K. Averaging over Uf yields 
now two terms, proportional to Cq-^+q2 and Cq-^^q^, re- 
spectively. While the first one is non-zero only for 
91+92 = 2 traversals of the link as required, the sec- 
ond one seems to spoil the proof. However, summing 
over qi at fixed 91+92, we find that such terms cancel in 

view of J2q Cq = C_2 + Cq + C2 = -1/2 + 1 + (-1/2) = 0. 

Having established that only paths visiting each link or 
2 times are to be considered, we can express, in analogy 
with 1^, the correlation functions (H) and (Q) in terms 
of sums over paths representing hulls in the percolation 
problem. In particular, we get for the products of Green 
functions entering (|^) 

(TrG(e', e; z)G{e, e'; z)) = (TrG(e', e; z-'^)G{e, e'; z'^)) 
= -2 5]P(e',e;iV)z2^, (8) 

N 

(TrG(e',e;z)G(e,e';z-i)) = -2 ^ Pi(e', e; TV)^^^, (9) 

N 

where P{e',e;N) and Pi{e' ,e; N) are probabilities that 
the edges e and e' belong to the same loop of the length N 
(resp. with the length N of the part corresponding to the 
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motion from e to e'). At z = 1 both (g) and (|^) reduce 
to —2 times the probabiUty P{e' , e) that e and e' belong 
to the same loop, which is identical (up to the sign) to 
the expression for the average conductance {g{e' , e)) at 
e = obtained in 
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FIG. 2. Scaling of the two-point conductance with dis- 
tance r between the contacts: average value (empty sym- 
bols), ((?), and typical value (filled symbols), gtyp — exp{lng), 
in systems of sizes L = 128(D) and L = 196(o). Also 
shown is scaling of the two-point Green function, (|Gp) 
and |G|?yp = exp(ln|Gp) (L = 128(A), L = 196(o)). The 
lines correspond to the r~^^^ (dotted) and r~^^'^ (dashed) 
power laws. Deviations from power-law scaling at large 
values of r are due to the finite system size. Inset: 
Fourier transforms of the one- and two-eigenfunction cor- 
relation functions, Si{r) = {\'ipia{e)i>ii3ie')\) (upper curves) 
and S2{r) = (7/)*„(e)7/)j-c(e)V'i/3(e')t/)*^(e')> (lower curves) for 
tij ~ (5 and L = 128 (solid), 256 (dashed), 384 (dotted). The 
dot-dashed line indicates a power law S{k) oc A;"^''* corre- 
sponding to S{r) oc r"^''^. 

The fractal dimension of the percolation hulls is 7/4 
10, implying that P and Pi scale as 

P(r, N), Pi{r, N) - N'^/'^r-^/^ , r < N^/^ (10) 

and fall off exponentially fast at r 3> N^/'^ (r is the dis- 
tance between e and e'). This yields for the correlation 
functions in (||) and (^) (which we abbreviate as {GrGr), 
{GaGa), {GrGa)) 



{GrGr) = {GaGa) ^ (GrGa) ^ r 



'1/2 



(11) 



in full agreement with the scaling argument of ^] 
and with our numerics. However, when we substitute 
®- ® ™ ®' these leading order terms cancel since 
^^[P(r, N) — Pi(r, N)] = 0. The result is non-zero due 
to the factors only, implying that relevant N are 
now N ^ 7~^, so that {{Gr — Ga){Gr — Ga)) scales 
differently compared to (pd]). 



I?(e',e;7) - ^ ^[P(r, iV) - Pi(r, A^)](I - g-^^^) 



N 



P{rn-')l~' - {i,rr^'\ r<i,. (12) 



Using now the definition 

-1/4 



p{e) ^ e^/'' ~ we find for r < ^, 



of V and the DOS scaling, 
(13) 



The same scaling behavior is obtained for the correlation 
function (|V'i^Q(e)'i/'j^(e')|) We thus conclude that 

rj =^ 1/4, consistent with our above numerical results. To 
shed more light on the difference in scaling between (GG) 
(or (g)), Eq. (pi]), and V, Eq. ( |l2|) , it is instructive to re- 
verse the logic and to ask how pi|) can be obtained from 
the wave function correlations (|12|) , (|l^) . It is straight- 
forward to express (GG) through V in the form of an 
integral over ei 2 with corresponding energy denomina- 
tors (dispersion relation). The integral is then dominated 
by £1,2 ~ £(?"), where e(r) is defined by ^^(^r) ~ 1^ (*-e- 
e(r) ~ This yields (GG) ~ V{r;e{r)) - r-^/^ 

(we have used (^ in the last step), in agreement with 
([ill). Therefore, (GG) (or (g)) is determined by wave 
functions with energies e(r), which transforms in 
( |T^ ) into an additional factor r~^/^. We will come back 
to this argument below to get an analogous relation for 
higher moments. 

In order to study the whole multifractal spectrum 
Aq, we return to numerical simulations. Our procedure 
based on the evaluation of ensemble averaged moments 
(IV'^'^I) and extrapolation to L ^ cx) was described in 
detail in |jl^. The results for arc shown in Fig. ||a. 
The obtained spectrum is parabolic with a high accu- 
racy. A parabolic spectrum is uniquely determined by rj, 
Aq ~ 7717(1 — q)/2; the above result rj = 1/4 thus implies 



A, ~ g(l - (7)/8. 



(14) 



We find, however, clear deviations from the parabolic law 
(p^), as shown in Fig. ^jb. One could ask whether these 
are not an artifact of uncontrollable finite-size corrections 
to scaling. We observe, however, an almost perfect scal- 
ing for all the quantities studied, yielding, in particular, 
?7 — 0.252 ± 0.002 in a very good agreement with the 
exact value 77 = 1/4. On the other hand, we find (in 
standard notations, aq — drq/dq) ao — 2 — 0.137 ± 0.003 
and 2-ai = 0.130 ± 0.003 (Fig. |c), while both these 
quantities would be equal to 77/2, should the parabolic 
law (|lj) be exact. This makes us believe that the spec- 
trum is only approximately parabolic, in contrast to ex- 
act parabolicity found for the QH transition |p^ . 

Finally, we turn to the statistics of two-point conduc- 
tances. Generalizing the above argument, we get 



([TrG(e',e)G(e,e')]') 



{e{r))L^HWll{e)i,l:^{e') 
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with an index Xg related to Ag and to the scahng dimen- 



sion Xp of the DOS (defined by p{e) 
(see also Jl6[), 

Xq = 2qxp + 2Ag. 
Using Xp — f/4 and Eq. (|lj) for A^, we find 



X„ 



q{3 - q)/4. 



) as follows 
(15) 

(16) 

The same scaling is expected to hold for the two-point 
conductance, {g'^) ~ r"'^'. However, since g is bounded 
from above, g <2, the exponent for {g'^) should be a non- 
decreasing function of q. Therefore, Eq.(16) will hold 
only for the moments (g*^) with q < qc, where qc — 3/2 
is the maximum of (|l6|), while for higher q the expo- 
nent saturates, Xq>q^ — Xq^ ~ 9/16 (these moments 
are determined by the probability to find g ^ 1). Equa- 
tion ( p^ ) corresponds to a normal distribution of In g (at 
r ^ 1) with the average (lnf;(r)) ~ — | Inr and the vari- 
ance var[ln g{r)] ~ 5 Inr. These arguments are fully con- 
firmed by the results of the numerical simulations shown 
in Fig. |. 
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FIG. 3. a) anomalous dimension A, extrapolated from en- 
sembles with system sizes L — 16 — 384. The dashed line 
(almost indistinguishable) is the parabola A, = 5(1 — g)/8; 
b) Aq/q{l — q) over q highlighting weak non-parabolicity; c) 
multifractal exponents a,(L) = -(|V'|^'' In |?/>p)/(|?/>p'') In L 
with q — (o), 1 (□), and 2 (o), and extrapolation to in- 
finite system size (dashed lines); ao ~ 2 — 0.137 ± 0.003, 
2 - ai = 0.130 ± 0.003, (2 - Q2)/3 = 0.125 ± 0.001. The 
arrow indicates the value ri/2 — 1/8 of these quantities for a 
parabolic spectrum. 

It remains an open question whether the multifractal 
exponents Ag, Xq can be computed by the conformal 
field theory methods |8|jl^Jl§^. Note that our results 
do not confirm the proposal of Ref . jl^] , where the value 
r] = 1/2 was obtained. Apparently, this indicates that the 
theory considered in |^ and obtained ||l^ from a partic- 
ular network model with fine-tuned couplings, does not 
belong to the SQH universality class. 



To summarize, we have studied, by combining numer- 
ical and analytical methods, the wave function statistics 
at the SQH transition. In particular, we have shown, 
using a mapping to classical percolation, that the index 
77 = — A2 [defined by Eq. (^] is equal to 1/4, at vari- 
ance with the r^^/^ scaling of the diffusion propagator 
n — {GrGa)- The multifractal spectra of wave functions 
(Ag) and two-point conductances {Xq) are given with a 
good accuracy by Eqs. ( p^ ) and (|l^) but show detectable 
deviations from parabolicity. 
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